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Abstract 

We study the static three-quark (3Q) potential in detail using SU(3) lattice QCD 
with 12^ X 24 at /3 = 5.7 and 16^ x 32 at /3 = 5.8,6.0 at the quenched level. For 
more than 200 patterns of the 3Q systems, we numerically derive 3Q ground-state 
QQ I potential V3Q from the 3Q Wilson loop with the smearing technique, which reduces 

f^ i excited-state contaminations. The lattice QCD data of V3Q are well reproduced 

^3 I within a few % deviation by a sum of a constant, the two-body Coulomb term and 

C— ^ . the three-body linear confinement term tragLniin, with Lmin the minimal value of 

^""^ ' total length of color flux tubes linking the three quarks. From the comparison with 

the Q-Q potential, we find a universality of the string tension as a^Q ~ ctqq and the 
one-gluon-exchange result for Coulomb coeflicients, ^3q ~ |^qq- 

Ph! 1 Introduction 

Strong interaction in hadrons or nuclei is fundamentally governed by quantum chromo- 
dynamics (QCD). In a spirit of the elementary particle physics, it would be desirable to 
^\f ' construct hadron physics and nuclear physics at quark level based on QCD. For instance, 

to bridge between QCD and hadron physics, the inter-quark potential[j|, ^, |3| is one of 
the most important quantities, because it is directly responsible for the hadron structure 
at the quark levelQ. However, it still remains as a difficult problem to derive the inter- 
quark potential analytically from QCD, because of the strong-coupling nature of QCD 
in the infrared region. 

Recently, the lattice QCD calculation has been adopted as a useful and reliable 
method for the nonperturbative analysis of QCD. For instance, quark-antiquark (QQ) 
potential, which is responsible for the meson properties, has been well studied using lat- 
tice QCD||l], |5|. In Fig.l, we show the QQ static potential UQQ(r) as the function of the 
distance r between quark and antiquark in SU(3) lattice QCD at the quenched level. 
Here, the QQ potential can be well reproduced by a sum of the Coulomb term due to 
the perturbative one-gluon-exchange (OGE) process, a linear confinement term and a 
constant ||l], ^, 

^QQ = -^ + ^QQ^ + Cqq, (1) 

where the string tension ctqq ~ 0.89 GeV/fm represents the confinement force. The 
linear potential at the long distance can be physically interpreted with the fiux-tube 
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Figure 1: The QQ static potential VQQ^r) obtained in SU(3) lattice QCD at the quenched 
level. 



picture or the string picture for hadronsH, 0, ^, P]. In this picture, the quark and 
the antiquark are linked with a one-dimensional flux-tube with the string tension ctqq. 
Hence, QQ potential in the long-range region is proportional to the distance r between 
quark and antiquark. This flux-tube picture (or the string picture) in the infrared region 
is supported by the Regge trajectory of hadrons, empirical analysis of heavy quarkonia 
data, the strong-coupling expansion of QCD0 and recent lattice QCD simulations. 

However, there is almost no reliable formula to describe the three-quark (3Q) potential 
VsQ directly based on QCD, besides the strong-coupling QCD[^, 0, although Vsq is 
directly responsible for the baryon properties[Q, ^, |^]. Up to now, the 3Q potential has 
been treated phenomenologically or hypothetically more than 20 years. We carry out the 
detailed study of the 3Q potential using SU(3) lattice QCD. 



2 Theoretical Consideration 

Theoretically, the 3Q potential[Q, |^ ^ is expected to be expressed by a sum of a constant, 
the two-body Coulomb term from the perturbative OGE process at the short distance 
and the three-body linear confinement term at the long distance, similarly in the QQ 
potential. Here, refiecting the SU(3)c gauge theory in QCD, the color fiux tube has a 
junction which combines three different colors, R, B and G, in a color-singlet manner. 
For most 3Q systems, the flux-tube energy is minimized with the presence of the junction, 
and therefore a Y-type flux tube is expected to be formed among the three quarks ||l|, Q, 0, 
P, 10|. We show in Fig. 2 the flux-tube configuration which minimizes the total fiux-tube 



length in the 3Q system. 

Thus, the static 3Q potential is considered to take a form as 
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+ CTSQ-^min + C'sQ, 



(2) 



where L^^^ denotes the miminal value of the total length of color fiux tubes linking three 
quarks. 




Figure 2: The flux-tube configuration of the 3Q system with the minimal value of the 
total flux-tube length L^in- 



Denoting by a, b and c the three side lengths of the 3Q triangle as shown in Fig. 2, 
Lmin is explicitly expressed as 111 



Lr, 



-(a^ + 6^ + c^) + -— ^'(a + b + c){-a + h + c){a -b + c){a + b-c) 



(3) 



when any angle of the 3Q triangle does not exceed 27r/3. For the case with one angle 
larger than 27r/3, Lmin is given as Lmin = (a + 6 + c) — max(a, 6, c). 



3 Static 3Q Potential in SU(3) Lattice QCD 

The static 3Q potential can be extracted from the 3Q Wilson loop [Q, ^, ^ in SU(3) 
lattice QCD calculations. As shown in Fig.^ the 3Q Wilson loop physically represents 
the creation of a gauge-invariant 3Q state at t = 0, the system with spatially fixed three 
quarks in < t < T and the annihilation of the 3Q state at t = T. 
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Figure 3: The 3Q Wilson loop in the Euclidean space-time. Tk{k = 1,2,3) denotes the 
path linking O and O' . 



In four-dimensional Euclidean space-time, the 3Q Wilson loop is a color current de- 
fined in a gauge-invariant manner as 



W^3Q = ^.eabcSa'h'c' uf U^'' uf , Uk = P expji^ / dx^^A^^x)] G SU(3) . 
3! JPfc 



(4) 



Here, P denotes the path-ordered product along the path denoted by Tj.{k = 1,2,3) 
linking O and O' in Fig.^. Similar to the derivation of the QQ potential from the Wilson 
loop, the 3Q potential Vsq can be obtained as Vsq = — limr^oo y ln(W3Q). 



Next, let us consider the physical state of the 3Q system. The ground state of the 
3Q system is expected to be composed by flux tubes rather than the strings|j|, P, and we 
denote by |g.s.;t) the 3Q ground state at t. On the other hand, there are many excited 
states of the 3Q system corresponding to the flux-vibrational modes[||], and we denote 
by |nthe.s.;t) the n-th excited 3Q state at t. In the 3Q Wilson loop, the gauge-invariant 
3Q state created at i = and annihilated at t = T can be expressed as 



|3Q;0) = co|g.s.;0) + ci|lst e.s.;0) +C2|2nd e.s.;0) + ..., 
|3Q;r) = co|g.s.;r) + ci|lst e.s.;r) + C2|2nd e.s.;r) + ... 

2 



(5) 
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where the coefficients are normalized as J2i=o IC' 
be consistent with the definition of W^q in Eq 
can be expressed as 



= 1. (This normalization is found to 
Then, the expectation value of W3Q 



(Wsq) = |co|2(g.s.;r|g.s.;0) + |ci|2(lste.s.;r|lste.s.;0) + ... 
= |co|^exp(-T/g.s.r) + |cipexp(-Vist e.s.T) + ... 



(6) 



with the ground-state potential V^.s. of the 3Q system and the n-th excited-state potential 

•"^n-th e.s. • 

As increasing T, the excited-state components drop faster than the ground-state com- 
ponent in (Wsq), however, the ground-state component |copexp(— V^.s.T) also decreases 
exponentially. Hence, we face a practical difficulty in extracting the numerical signal. To 
avoid this difficulty, we adopt the smearing technique |]I|, P, |ll| which enhances the ground- 
state overlap and removes the excited-state contamination efficiently. (There were a few 
lattice studies on the 3Q potential more than 13 years ago[]T2|, p^ , however, their results 
seem unreliable and misleading due to the fatal large excited-state contaminations. The 
smearing technique was mainly developed after their works.) 

4 Smearing Technique to Remove Excited Modes 





Figure 4: The schematic explanation of the smearing for the link- variables. 



The smearing for the link-variable is expressed as the iterative replacement of the 
spatial link- variables by a obscured link- variables Ui{s) £ SU(3) which maximizes 



Re tr i U}(s 



a 



Ui{s) + J2{Uj{s)Ui{s + j)C/j(s + i) + Uj{s - j)Ui{s - j)Uj{s + i - j)} 



with a real smearing parameter a. This can be visualized as in Fig.^. The obscured 
spatial line composed by the obscured link-variables physically corresponds to the flux 
tube with a finite cylindrical radius. The smearing parameter a and the iteration number 



(7) 



of the smearing are suitably determined so as to maximize the ground-state overlap of the 
3Q system at i = and i = T in the 3Q Wilson loop. In fact, after a suitable smearing, 



we expect saturation of the ground-state overlap as |co| 
the excited-state contaminations as |cjp ~ (i = 1,2,..) 



~ 1, i.e. strong reduction of 
J, g, 11|. To see this in actual 
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Figure 5: The ground-state overlap Co = {W3q{T))^+'^ / {W3q{T + 1))'^ in SU(3) lattice 
QCD at /? = 5.7. By the smearing, the ground-state overlap Cq is largely enhanced from 
the lower data to the upper data as 0.8 < Cq < 1 for each 3Q system. 



lattice QCD calculations, we examine in Fig.|5| the ground-state overlap 

Co^{Wsci{T)f+'/{Wsci{T+l)f 



(8) 



corresponding to |coP[T|, ^], and we observe a large enhancement of the ground-state 
overlap as 0.8 < Cq < 1 by the smearing. To summarize here, owing to the smearing, we 
can set up the ground-state dominant 3Q state at t = and t = T in the 3Q Wilson loop, 
and therefore we enjoy accurate measurements for the 3Q potential without suffering from 
excited modes like the flux-tube vibration. 



5 The Lattice QCD Results for The 3Q Potential 



We study the 3Q potential V3Q in detail by investigating more than 200 patterns of the 
3Q systems using SU(3) lattice QCD with 12^ x 24 at /3 = 5.7 and 16^ x 32 at /? = 5.8, 6.0 
at the quenched level. Here, the lattice spacing a is determined so as to reproduce the 
QQ string tension CTqq = 0.89GeV/fm at each p. We show in Fig.^ the 3Q potential 
V3Q in SU(3) lattice QCD at 13=5.7 and /3=5.8 as the function of Lmim the minimal 
value of the total flux-tube length. At the large distance, where the Coulomb potential is 
negligible, we observe the linearity on Lmjn, which is consistent with the expected form in 
Eq.(|2|). At the short distance, perturbative QCD would be applicable, and therefore the 
two-body Coulomb potential in V3Q can be estimated as Vcoui 
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the Coulomb coefficient ^qq in the lattice QCD result of the QQ potential. (Here, the 
color factor ^ appears in perturbative QCD.) To single out the confinement contribution, 
we subtract this perturbative Coulomb contribution Vcoui from the 3Q potential V3Q. 



The results are shown in Fig.0. The clear linearity on Lmm in the whole region observed 
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Figure 6: The 3Q potential V3Q as the function of i^mim the minimal value of the total 
flux-tube length, in SU(3) lattice QCD at /? = 5.7 (left) and at /? = 5.8 (right). 
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Figure 7: The SU(3) lattice QCD result for V3Q — Vcoui as the function of Lmin, the 
minimal value of the total flux-tube length, at /3 = 5.7 (left) and at /3 = 5.8 (right). 



in Fig.0 means that the 3Q potential V3Q can be well reproduced by a sum of the the 
perturbative Coulomb contribution as Vcoui and the linear confinement term proportional 
to Lmin in Eq.dl). 

Now, we consider the best fitting of the 3Q potential data with the form in Eq.(g). 
With three parameters A^q, ctsq and C3Q, the 3Q potential data can be well reproduced 
within the accuracy better than a few %. We show the best fitting parameters of the 3Q 
potential V3Q and the QQ potential Vqq in Table [l| at various /?. 

First, we focus on the properties at the short distance, where the nonperturbative 
confinement force is almost negligible. We find A^q ~ 2^QQ ^^^ ^^^ Coulomb coefficient, 
and find the perturbative color factor ^ originating from the OGE contribution. 

Second, let us consider the constant term. The constant terms C3Q and Cqq behave 
as -, and diverge in the continuum limit a ^ 0, so that the constant term is not physical 
quantity. In the limit of r ^ in the QQ system, the Wilson loop becomes trivial 
as WQQ{r = 0) = 1 and then Vqq(?^ = 0) = holds, while one gets VQQ{r ~ 0) ~ 

— -^-I-Cqq with < a; < 1 [|l], 1^ in the lattice regularization. Then, we find Cqq ~ -^. 



From the similar argument on the diquark hmit, we find also C3Q — Cqq ~ -^ [||, |^. In 
this way, the constant term reflects the Coulomb singularity, and is kept to be finite in 
lattice QCD, because of the lattice reguralization[|^, |3[. Physically, these constant terms 
C3Q and Cqq represent the pairing energy at the limit when the quarks gathers at one 
spatial point. From the above two relations and A^q ~ ^^qq) we obtain C3Q ~ I^qq, 
which seems to hold as shown in Table 1. 

Finally, we consider the infrared nonperturbative properties of the inter-quark po- 
tential. The universality of the string tension is observed as cjsq ~ '^qq- ^^ fact, the 
strength of the confinement force is universal for the QQ system and the 3Q system, 
which supports the color flux-tube picture for hadrons y, |8|, [3, |§, P^] - 



6 Summary and Concluding Remarks 

To bridge the gap between QCD and hadron physics, we have studied the 3Q potential 
responsible for the baryon properties. We have calculated the ground-state 3Q potential 
for more than 200 patterns of the 3Q systems using SU(3) lattice QCD at /3 = 5.7, 5.8 and 
6.0 at the quenched level. In this calculation, we have adopted the smearing technique 
which enhances the ground-state overlap and reduces the excited-state contaminations 
like the flux vibrational modes. Within the accuracy better than a few %, the lattice 
QCD data for the 3Q potential have been well reproduced by a sum of a constant, the 
two-body Coulomb term, and the three-body linear confinement term proportional to 
-^miii) the minimal value of the total flux-tube length. We have observed the perturbative 



OGE result as A 



3Q 



— 2^C 



j, and the universality of the string tension as (T3Q ~ o"qq- 
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a (GeV/fm) 


A 


C (lattice unit) 


3Q (/3 = 5.7) 
QQ (/? = 5.7) 


0.832(15) 
0.890(25) 


0.1331( 66) 
0.2793(116) 


0.9182(213) 
0.6203(161) 


3Q (/3 = 5.8) 
QQ (/? = 5.8) 


0.818( 6) 
0.890(21) 


0.1304( 17) 
0.2580(159) 


0.9326( 53) 
0.6081(182) 


3Q {13 = 6.0) 
QQ {f3 = 6.0) 


0.811( 7) 
0.890(12) 


0.1363( 11) 
0.2768( 24) 


0.9590( 35) 
0.6374( 30) 



Table 1: The best fitting parameters for the 3Q potential V3Q at /3=5.7, 5.8 and 6.0. 



Table 2. Partial data for the 3Q potential in SU(3) lattice QCD at (3—5.8. {i,j, k) labels the 
3Q system where the three quarks put on (i, 0, 0), (0, j, 0) and (0, 0, k) in R'^ in the lattice unit. 



J, 



T/latt 
^3Q 



0.9144 
1.0657 
1.1962 
1.3049 
1.4134 
1.5339 
1.1866 
1.3062 
1.4153 
1.5259 
1.6360 
1.4178 
1.5216 
1.6276 
1.7296 
1.6291 
1.7190 
1.8203 
1.2793 
1.3815 
1.4904 
1.4751 
1.6814 



^3Q 



0.9008 
1.0582 
1.1883 
1.3054 
1.4162 
1.5236 
1.1877 
1.3076 
1.4203 
1.5287 
1.6347 
1.4209 
1.5296 
1.6357 
1.7401 
1.6356 
1.7396 
1.8425 
1.2838 
1.3903 
1.4969 
1.4880 
1.6912 



T/latt 
^3Q 



'^3Q 



0.0136 
0.0075 
0.0079 
-0.0005 
-0.0029 
0.0103 
-0.0011 
-0.0014 
-0.0050 
-0.0029 
0.0013 
-0.0031 
-0.0080 
-0.0081 
-0.0105 
-0.0064 
-0.0206 
-0.0222 
-0.0046 
-0.0089 
-0.0064 
-0.0129 
-0.0098 



{hJ,k) 


T/latt 
'^3Q 


^3Q 


T/latt T/tit 
^3Q ~ ^^3Q 


(4,5,6) 


2.1263 


2.1094 


0.0169 


(4,6 


6) 


2.2264 


2.1979 


0.0285 


(5,6 


6) 


2.3105 


2.2733 


0.0372 


(0,1 


1) 


0.7701 


0.7727 


-0.0027 


(0,1 


2) 


0.9659 


0.9694 


-0.0035 


(0,1 


3) 


1.1107 


1.1071 


0.0037 


(0,1 


4) 


1.2323 


1.2268 


0.0055 


(0,1 


5) 


1.3339 


1.3388 


-0.0049 


(0,1 


6) 


1.4566 


1.4469 


0.0097 


(0,1 


7) 


1.5730 


1.5527 


0.0203 


(0,1 


8) 


1.6779 


1.6570 


0.0209 


(0,2 


2) 


1.1375 


1.1422 


-0.0047 


(0,2 


3) 


1.2660 


1.2711 


-0.0050 


(0,2 


4) 


1.3815 


1.3871 


-0.0056 


(0,2 


5) 


1.4928 


1.4972 


-0.0043 


(0,2 


6) 


1.6069 


1.6041 


0.0028 


(0,2 


8) 


1.8285 


1.8128 


0.0157 


(0,3 


3) 


1.3893 


1.3940 


-0.0047 


(0,3 


41 


1.5007 


1.5066 


-0.0059 


(0,3 


5) 


1.6135 


1.6145 


-0.0010 


(0,4 


4) 


1.6079 


1.6165 


-0.0086 


(0,4 


5) 


1.7085 


1.7226 


-0.0141 


(0,4 


6) 


1.8263 


1.8267 


-0.0004 
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